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ABSTRACT 

We report the discovery of a short-duration microlensing candidate in the 
northern field of the POINT- AGAPE pixel lensing survey towards M31. The 
full-width half-maximum timescale is very short, ^1/2 = 1.8 days. Almost 
certainly, the source star has been identified on Hubble Space Telescope archival 
images, allowing us to infer an Einstein crossing time of = 10.4 days, 
a maximum magnification of ^max ~ 18, and a lens-source proper motion 
/Xrei > 0.3 //as/day. The event hes projected at 8' from the center of M31, which 
is beyond the bulk of the stellar lens population. The lens is likely to reside in 
one of three locations. It may be a star in the M31 disk, or a massive compact 
halo object (Macho) in either M31 or the Milky Way. The most probable mass 
is 0.06 M0 for an M31 Macho, 0.02 M© for a Milky Way Macho and 0.2 M© for 
an M31 stellar lens. Whilst the stellar interpretation is plausible, the Macho 
interpretation is the most probable for halo fractions above 20%. 

Subject headings: Galaxy: halo - M31: halo - lensing - dark matter 
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1. Introduction 

Following the suggestion of Paczyhski (1986), several groups searched for dark matter 
in the form of massive compact halo objects (Machos), using gravitational microlensing of 
background stars in the Magellanic Clouds (see Alcock et al. 2000; Lasserre et al. 2000; 
and references therein). After monitoring about 10'' stars for several years, the results are 
consistent with a Macho halo mass fraction of 20%, though with considerable uncertainty. 

The possibility of detecting microlensing events in M31 was independently suggested 
by Crotts (1992) and Baillon et al. (1993). The advantage of targeting a large external 
galaxy is that the number of stars that can act as possible sources is enormous. Moreover, 
the high inclination of M31's disk causes an asymmetry in the observed rate of microlensing 
by lenses in a spheroidal halo. Although this gives an unambiguous signature of the halo 
lenses (|Crotts 1992|) , the difficulty is that the sources are resolved only while they are lensed 



(and then only if the magnification is substantial). Nonetheless, pilot campaigns by the 
Columbia- VATT group ( |Crotts fc Tomaney 1997| ) and by the AGAPE collaboration ( |Ansari 



et al. 19971 ; [Ansari et al. 19991) established the feasibility of such observations and identified 
some candidate microlensing events. In particular, the AGAPE collaboration ([Ansari et al 
1999) reported a short duration candidate, AGAPE Zl, in the bulge of M31. 



The POINT-AGAPE collaboration^ employs the Wide Field Camera (WFC) on the 
2.5 m Isaac Newton Telescope (INT) to carry out a pixel-lensing survey towards M31, 
monitoring two fields of 0.3 deg^ each, located North and South of the M31 center. The 
survey has the potential to map the global distribution of the microlensing events in M31 
and to determine any large-scale gradient. POINT-AGAPE is a long-term program, as at 
least three years of data will be required to make a convincing identification of a gradient 



^ see [http : / /www . point-agape . org 
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(Kerins et al. 2001). Also, it is usually necessary with pixel lensing to establish a long 
baseline to distinguish microlensing events from variable stars. In this Letter, we report on 
a particularly interesting and convincing early candidate, PA-99-N1, so named because it is 
the first event to be announced by POINT-AGAPE that peaks in the northern field in 1999. 



2. Observations and Data Analysis 

We restrict analysis to a (11' x 22') field centered 3' west and 12' north of the center of 
M31. The observations are spread over 36 epochs between August and December 1999. The 
exposures are in two bands: 36 epochs in Sloan r' and 26 epochs in Sloan g'. The exposure 
time is typically between 5 and 10 minutes per night. The sampling averages to one epoch 
every three nights, though the observations are strongly clustered because the WFC is not 
always mounted. The data reduction is described in detail elsewhere ( |Ansari et al. 1997| ; 
Strauman et al. 19981 ; |Le Du 20001) . 



After bias subtraction and fiat-fielding, each image is geometrically and photometrically 
aligned relative to a reference image (14 August 1999) which was chosen because it has 
a long exposure time, typical seeing (l."6) and little contamination from the Moon. The 
lightcurves are computed by summing the flux in 7-pixel (2."3) square "superpixels" and 
then removing the correlation with seeing variation. This size is set by the worst seeing 
~ 2."1. 

We use a simple set of candidate selection criteria designed to isolate high signal-to-noise 
events. Detection of events is made in the r' band, which has better sampling and lower 
sky background variability. For an event to be detected, it must induce one and only one 
significant bump on the lightcurve. A bump is defined by at least 4 consecutive data points 
rising above the background by a minimum of 4a. Its significance is quantified by the 
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probability P of obtaining data points, each of which has at least the signal to noise ratio 
of the measured points (evaluated assuming Gaussian errors). To microlensing 
candidate, we demand the lightcurve has only one bump with — InP > 100 and no other 
bump with — InP > 20. A similar selection procedure was used by the AGAPE Pic du Midi 
program (Ansari et al. 1997). 

The background is the minimum of the running average of 7 consecutive data 
points. Flagged events are fitted to a high magnification degenerate microlensing 
curve (Gould 1996) simultaneously in the r' and g' bands, whose 6 parameters are 
(ti/2,to,-Fbase,r',-Fbase,5', AFpeak,r', AFpeak,g')- Here, ti/2 is the full-width at half-maximum, to 
is the time of the peak, -Fbase is the baseline flux and AFpeak is the flux difference between 
baseline and maximum. We require a P^r degree of freedom be less than 3 for this fit. 
We also require at least minimal detection in g' , i.e. AFpeak.g' > 1 ADUs^^. To insure a 
high probability of detecting microlensing rather than other forms of stellar variability, 
given the clustered sampling and relatively short duration of the observations, we demand 
ti/2 < 8 days and a peak flux AFpeak,r' > 10ADUs^\ corresponding to -Rpoak < 21.5. We 
find one candidate which has ti/2 ~ 1.8 day and AFpeak,r' ~ 17ADUs~^. 



3. The Microlensing Candidate 

Figure shows the lightcurves in r' and g' of this candidate together with the 
non-degenerate fit derived below. The g' data with comparatively large error bars were 
taken on nights with high Moon background. Using the Aladin Sky Atla^, we find that 
PA-99-N1 has J2000 position: a = 00^42°^51!42, 6 = +41°23'53'.'7. That is, it lies projected 
on the near disk, 7'52" from the center of M31. 



see [http : / /aladin . u-strasbg . f r/aladin . gml 
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There are some straightforward tests to see if PA-99-N1 is compatible with microlensing. 
First, there are no comparable "bumps" in the remainder of the lightcurve shown in Figure 
ID as might be expected for many classes of variable stars. (There is, however, a much 
smaller but somewhat disturbing bump in December, near day 130. We investigate this 
bump in § 4.1.) There also are no comparable bumps at this field position in our data taken 
in 1998 and 1999 with the 1.3m McGraw-Hill telescope at the MDM observatory. Second, 
microlensing events are achromatic, so the ratio of flux change in different bands should be 
constant in time (e.g., Ansari et al. 1997). This requires 

AFg/ (t) _ Fg' (t) — Fbase,g' 
AFr' (t) Fr' (t) — -Fbase,r' 

which does indeed hold for PA-99-N1 (see Fig. |I]c). 



constant, 



Using DAOPHOT ([Stetson 1987|) on the images taken nearest maximum magnification, 
we find R = 20.80±0.13 and V = 22.00±0.17, i.e., V-R = 1.20±0.22 and My = -2.8±0.3, 
assuming (m — M)o,m3i = 24.43 and estimated total extinction Ay = 0.4 ± 0.2. The 
transformation from instrumental (r', g') to Johnson {V, R) is based on 31 standard stars 
that lie in the same field as PA-99-N1 ( |Magnier et al. 1992| ; [Magnier et al. 1993 ; Haiman et 
al. 19941) . 



The PA-99-N1 position lies within a series of five Hubble Space Telescope (HST) 
WFPC2 archival images taken in July 1996, three with F814W (~ /) and two with F606W 
(~ V). We use the relations of Zheng et al. (2001) to transform from these filters to 
Johnson-Cousins V and /. Since the event is very red, V — R ^ 1.2, the source must lie 
either high up on the giant branch where essentially all stars are resolved by HST, or must 
be on the main-sequence and so magnified by Amax ^ 10'' (see Fig. ||). The latter possibility 
is extremely unlikely a priori. Hence, a firm prediction of the microlensing interpretation 
is that there should be a resolved star in the HST image with the same color and same 
position as the event. 
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To measure the spatial position of the event within the INT image we take the 
difference of the image at maximum and another image at basehne that has similar seeing. 
The difference image shows a well-defined and isolated difference star whose position we 
measure using IRAF IMEXAMIN. Variation of the IMEXAMIN parameters changes the 
result by ~ 0.2 pixels (0."07), which we adopt as the error in the measurement. The spatial 
transformation between the INT and HST fields is derived by comparing the positions of 
the eight stars that are both resolved in the INT image and unsaturated in the HST image. 
The uncertainty in this transformation is much smaller than that of the INT position 
measurement. Figure |^ shows part of the HST field with the la and Sex positional error 
circle centered on the event. Close to this circle there is a resolved star for which we find 
V = 24.51 ± 0.12 and I = 22.41 ± 0.10, implying V -I = 2.10 ± 0.16. There are no other 
resolved stars within the 3 cr circle. For typical stellar populations this V — I is compatible 
with the V — R = 1.20 ± 0.22 measured for the event (Demarque et al. 1996; Yi, Demarque, 
& Oemler 1996). The prior probability to find such a star so close to the predicted position 
is only 3%. We conclude that the resolved HST object is almost certainly the source star of 
the event. 

4. Lightcurve Interpretation 

4.1. PA-99-N1 as a Variable Star 

Before analyzing PA-99-N1 as a microlensing event, we need to ask whether it 
is consistent with an interpretation as a variable star? The color and magnitude are 
compatible with a MO Mira at maximum. However, such Miras have periods greater than 
200 days and their flux at maximum does not vary as rapidly as the event (Allen 1999). 
PA-99-N1 is unlikely to be a dwarf nova. It is too bright at maximum to be in M31. If 
it is a Galactic dwarf nova, it must be either unusually faint or improbably far from us 
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(> 10 kpc). It is also unlikely to be a nova. The brightness decrease of PA-99-N1 just after 
maximum is very rapid with a rate of decline of about 0.7 magnitudes/ day. Capaccioli et 
al. (1989) studied the relation between the rate of decline and the magnitude at maximum 
for novae in M31 and found that the brighter the nova, the faster the decrease. The decline 
of PA-99-N1 would indicate a nova as bright as V = 16, which would imply Ay = 6, in 
stark contrast to the small extinction seen in Figure |^. To conclude, there is no type of 
stellar variability known to us that could generate the lightcurve of PA-99-N1. 

However, if the December bump near day 130 in Figure |l] were due to the event source, 
one would nevertheless have to conclude that the source was probably a variable. We 
therefore investigate this bump closely. For each of the nine December images, we find a 
corresponding image at baseline with comparable seeing, and subtract the two. Summing 
these difference images, we find a very clear stellar profile located l."l ± 0."1 south of the 
source. Hence, this bump is due to genuine stellar variability, which was close enough to 
contaminate our 2."3 superpixels, but is clearly distinguishable from the event source. We 
eliminate these contaminated December points from the analysis. 

4.2. PA-99-N1 as a Microlensing Event 

Making use of the identification of the HST star as the source of the microlensing 
event, we fit the data to a full (non-degenerate) microlensing curve and so evaluate the 
Einstein crossing time, = 10.4 days. 

By permitting the measurement of ^e, 



the identification of the source reduces but does not eliminate the degeneracy in the lens 
parameters. The mass M remains entangled with vTrei and firei, the lens-source relative 
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parallax and proper motion. If the lens transits the source, one can measure /irei and so 
further break the parameter degeneracy (Gould 1994; Alcock et al. 1997). We are not able 
to measure fi^ei for PA-99-N1, but are able to strongly constrain it from the lack of finite 
source effects that would be induced by a transit. 

Let p* = 6^:/6e, where 0* is the angular radius of the source. We find that if we fit 
the lightcurve with fixed at any value p* < 0.1, then the other parameters all assume 
the same values as in Table 1. However, as p^, is increased further, t-^ declines, so that 
the parameter combination p^t-^ saturates at 1.0 day, and then for p^, ;^ 0.12, rises 
dramatically. From this we derive the constraint p*tE < 1.0 day, implying, 

A^rel = = > • (3) 

tE pJe 1.0 day 

Using the empirical surface-brightness/color relation of van Belle (1999), the color-color 
transformation of Bessell & Brett (1988), and an estimated total extinction Ay = 0.4 ± 0.2, 
we find 6'^ = 0.30 ± 0.04 pas, where the uncertainty is dominated by the 11% intrinsic 
scatter about the van Belle (1999) relation. The constraint (§) then yields 



.„, = ,„,A> 400 Urns- (^) (- 



(4) 

.3 pas ) \^770 kpc j ' 

where frei is the transverse speed of the lens relative to the observer-source line of sight, and 
Di is the distance to the lens. While the constraint (^) is unimportant for Galactic lenses 
(where Di ~ 10 — 30 kpc), it strongly limits the allowed range of prei for M31 lenses (since 
high values are exponentially suppressed). Hence prci is crudely measured, which implies, 
via equation (0), a constrained relation between M and vTrei. Thus, when the proper-motion 
constraint @) is incorporated into the Monte Carlo (see below) it indirectly constrains 
the mass to a much narrower range than would be allowed without it. Note that our 
finite source fits incorporate a linear limb-darkening parameter of 0.78 in R band, which is 
appropriate for a V — / = 2.1 (MO-Ml, 3630 K) star (Manduca, Bell & Gustafsson 1977; 
Johnson 1966). 



5. Discussion 



The location of PA-99-N1, at nearly 8' from the center of M31, is interesting because 
the majority of stellar lenses are expected to reside within 5' of the M31 center (e.g., Kerins 
et al. 2001). The significance of the candidate is assessed by performing Monte Carlo 
simulations of events with the same source magnitude, Einstein timescale and projected 
position as PA-99-N1, as well as taking account of the lower limit on /irei in equation (^). 
We use the actual sampling and exposure times for the 1999 season. We model the halos of 
both galaxies with cored near-isothermal spheres, taking the mass of M31 as twice the mass 
of the Galaxy and assuming a core radius of 5 kpc for both galaxies. The M31 bulge follows 
Kent's (1989) axisymmetric model while the disk has a sech-squared profile (see Kerins et 
al. 2001 for details). We model the disk stellar lens masses based on the Galactic disk mass 
function (MF) of Gould, Bahcall, & Flynn (1997) corrected for binaries and extended down 
to 0.01 M0, and the bulge stellar masses using the Zoccali et al. (2000) Galactic bulge MF, 
similarly corrected and extended. 

The thick solid line in Figure ^ (a) shows the relative probability P that PA-99-N1 
is due to a Macho of mass M, assuming full Macho halos. The curve is normalized with 
respect to the probability that the lens is a star, P^. The thick dashed (dotted) curve 
gives the contribution to P{M) from M31 (Milky Way) Machos. In addition to this, we 
have plotted the M31 distribution before applying the proper motion cut of equation (|) 
(thin dashed line). The black horizontal line ('d/b') gives the fractional contribution of 
the disk lens and bulge source configuration to P^,, whilst the gray line ('b/b') shows the 
contribution of bulge-bulge lensing. The width of these lines reflects the dispersion in the 
logarithm of the stellar lens mass. 

Clearly, if the halos are full of Machos, the most probable interpretation is that the lens 
is a Macho with mass ~ 0.03 Mq. This is about five times more likely than PA-99-N1 being 
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due to a stellar lens. The lens is about equally likely to lie in the Milky Way or M31 halo, 
which is quite unexpected given that the line of sight passes only 2 kpc from the center of 
M31, but 8 kpc from the Galactic center. Normally this would cause the M31 probability 
to be much higher (as shown by thin dashed line). However, the fi^ei constraint (H) severely 
suppresses the M31 distribution at low masses, while leaving the Milky Way distribution 
virtually unaffected. This is because the average Milky Way Macho has a relative proper 
motion of ~ 10 /ias/day, whilst for the typical M31 Macho it is ~ 0.2 /las/day. The 
constraint also has the effect of displacing the peak of the M31 distribution towards higher 
mass. 

Assuming a logarithmic prior in the Macho mass gives M = GAtfg X 10-2 Mo for an 
M31 Macho and 2.ltli x lO^^ M© for a Milky Way Macho. The M31 and Milky Way 
distributions are both quite broad. As they are about equal in amplitude and are displaced 
from one another by nearly a decade in mass, their combined distribution is even broader, 
with a FWHM extending from the planetary to the stellar mass range. 

If the lens is a star, then the overall mass probability distribution is given by the 
thick solid line in Figure ^ (b). The thick dashed line ('d/b') shows the contribution from 
disk-bulge lensing whilst the thick dotted line ('b/b') shows the distribution for bulge-bulge 
lensing. The contribution of lensing involving disk sources is negligible. In the absence 
of the relative proper motion cut, the disk and bulge distributions would be given by the 
thin dashed and dotted lines respectively. Including these cuts results in a most probable 
mass around 0.2 Mq, making a rather plausible hypothesis that the lens is a low mass star. 
Logarithmic averaging of the solid line gives M = 0.27t|]:?^ Mq, i.e., only a factor 1.8 (Icr) 
uncertainty. 

We see that, if the lens is an M31 star, it most likely lies in the disk. This is a direct 
consequence of the fact that the source hes at the edge of the bulge, while the line of sight 
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passes through the near side of the disk. Comparing Figures |^a and ^b, it is clear that 
the mass of the lens is much better constrained if it lies in the M31 disk than the M31 
halo. This is because for M31 lenses, M oc (/irei^E-Dis)^ where D\s is the distance from 
the lens to the source (see eq. 0). For both populations, the product /irei^E is reasonably 
well constrained, but for halo lenses, D\s can take on a very broad range of values, while 
for disk lenses, the geometry implies Di^ = 4.0 ± 1.8 kpc. Since is measured and /irei is 
constrained by equation (^), M is also constrained. This means that we only detect events 
at the location and timescale of PA-99-N1 if the stellar mass M > 0.1 M0(4 kpc/Dis). 

The relative rate of Machos to stars is subject to a number of modeling uncertainties. 
First, we have assumed that the M31 halo is twice as massive as that of the Milky Way. 
The most recent mass estimates of the M31 halo derived from the kinematics of the satellite 
galaxies suggest that it may be only roughly as massive as the Milky Way (Evans et al. 
2000). Second, it is sensitive to the choice of the core radius for M31's halo, with larger 
core radii giving a reduced M31 Macho rate at the location of PA-99-N1. Third, we have 
assumed that the M31 bulge is axisymmetric, but the twisting of the optical isophotes 
(e.g., Walterbos & Kennicutt 1987) is evidence for the presence of a bar. Inspection of 
their Figure 5 shows that the twisting is away from the location of PA-99-N1. We therefore 
expect a bar model with the same overall mass as our axisymmetric model to have a lower 
surface density at the position of the event and thus a lower stellar lensing rate there. 



6. Conclusions 

We have reported the discovery of PA-99-N1, a high signal-to- noise, short duration 
event which is consistent with the microlensing hypothesis. Almost certainly, the source 
star has been identified on archival HST frames, from which the Einstein crossing time of 
10.4 days has been determined. We have argued - from the stability of the lightcurve, the 
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achromaticity of the flux excess, the excellence of the flt, and the consistency of the color 
of the event at maximum with the color of the HST source - that by far the most natural 
explanation of this event is microlensing. The lens is most likely to be a Macho if the halo 
fraction is above 20%. However, it is also plausible that the lens is a disk star with a mass 
of ~ 0.2 Mq. 
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Fig. 1. — Panels (a) and (b) show the flux in r' and g' against time in days. Panel (c) is a 
zoom centred on the event that shows the variation of the ratio of the flux change in the two 
passbands AFg/ / AFr> with time. The vertical lines are centred on to and are separated by 
0.9 days, i.e., half the full width at half-maximum. The days correspond to J — 2451392.5 
where J is the Julian date. 
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Fig. 2. — Superposition of the Hipparcos (points) and HST (shaded area) color-magnitude 
diagrams. The Hipparcos stars have been moved to the position of M31 [assuming 
(m — M)o = 24.43 and Galactic foreground extinction Ay = 0.24] . The intersection of 
the horizontal and vertical lines shows the position of the HST object. 




Fig. 3. — An image of the HST field showing the 1 cr and 3 a error circles around the position 
of the event. There is a resolved source close to the la circle. 
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Fig. 4. — Panel (a) shows, for the case where PA-99-N1 is a Macho, the probabihty P that 
the lens has mass M (thick solid line), normalized to the probability that PA-99-N1 is a 
stellar lens. The contribution to this curve from M31 (Milky Way) Machos is shown by the 
thick dashed (dotted) line. The thin dashed line shows P(M) for M31 Machos before the 
relative proper motion cut (eq. The horizontal bars show the contribution to from 
bulge-bulge lensing ('b/b') and disk-bulge lensing ('d/b'), whilst their width indicates the 
logarithmic dispersion in the stellar lens mass. Panel (b) shows, for the case where PA-99-N1 
is a stellar lens, the overall stellar probability density as a function of lens mass (solid line), 
including the contributions from bulge-bulge lensing (thick dotted line) and from disk-bulge 
lensing (thick dashed line). The thin dotted and dashed lines show the relative distributions 
of bulge and disk lenses before the proper motion cut. 
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Fit Parameters 


N 




to 


Uo 




Fs,r' 


Fs,g' 


Fb,r' 


Fb,g' 






(days) 




(days) 


(ADU/s) ( 


;adi:/s) 


(ADU/s) 


(ADU/s) 


97 


154 


13.87 


0.056 


10.4 


1.02 


0.28 


397.65 


209.40 




± 


0.04 


0.009 


1.5 


0.17 


0.05 


0.18 


0.10 



Table 1: Fit parameters for PA-99-N1. N: number of points, to- time of the peak, uq: impact 
parameter, Einstein timescale, Fg and F^: source and background fluxes in two bands. 
The maximum magnification is Amax ~ Mq ^ ~ 18 
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